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Funding Source 

• This work was funded by the NASA MSFC 
Space Launch System Advanced Development 
Project under Grant NNM13AA08G with Melinda 
Nettles as the Program Director 

• The purpose of this project is to develop new 
and improved two-phase flow boiling heat 
transfer correlations for the Generalized Fluid 
System Simulation Program (GFSSP) 
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Background 

In many cases, cryogenic systems 
require single-phase liquid transfer 
to a destination. 

For the cryogen to remain a liquid, 
the hardware must be chilled down 
to cryogenic temperatures. 

The amount of hardware to chill 
down can be quite large, such as a 
feed line to a propellant storage 
tank. 
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Background 

• Cryogen is fed through the hardware to chill it down. An initial portion of 
cryogen absorbs the heat from the pipe and boils off into vapor. 

• This quantity of cryogen is sacrificed to allow single-phase liquid transfer. 

• System design requires accurate knowledge of the necessary amount of 
sacrificed cryogen. Minimization of the sacrificed cryogen is desirable. 

• The amount depends on the q” from pipe to fluid and the FR of the fluid. 

• The robust, accurate prediction of the two-phase flow boiling heat flux 
continues to be a challenge 

n 

• h = = /(m, P, fluid properties, etc . ) 

T w -Tf 
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Boiling Regimes 


3 different boiling regimes 
during chilldown 

1. Film boiling : In this 
regime T w » T sat . Any 
liquid approaching the 
wall is vaporized 
before touching the 
wall. 
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Boiling Regimes 


3 different boiling regimes 
during chilldown 

2. Transition Boiling : T w 
drops below the 
rewetting 
temperature, T wet . 
Liquid is able to 
intermittently touch 
the wall. 
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Boiling Regimes 

3 different boiling regimes 
during chilldown 

3. Nucleate Boiling : T w drops 
below the temperature 
corresponding to the critical 
heat flux (CHF). This marks 
the start of full nucleate 
boiling. Liquid remains in 
contact with the wall. Heat is 
transferred by vapor bubbles 
generated in surface cavities 
and swept away from the 
surface. 
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Goals of this work 

• Determine accurate correlations for the film, transition, 
and nucleate boiling heat transfer coefficients 

• Determine accurate correlations for T wet , CHF, and 
T onb . These help determine the current boiling regime 
in a chilldown simulation. 

• Develop a 1-D numerical simulation of the chilldown 
of a pipe which uses these correlations to predict the 
wall heat flux and compare to the results to recent LN 2 
chilldown experiments 


8 


a& JUJk& = 

Shaping the Future of Aerospace 

Experiment 

• Wall temperature and wall heat flux measurements 
were obtained from experiments in which a thin, 
vertically-aligned tube was chilled down with the 
downward flow of LN 2 . 

• 54 tests were completed. 

• G = 61.2 to 1150 kg/m 2 s 

• P = 174 to 817 kPa 

• Subcooled inlet temp from 0 to 14 K 
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Experiment Schematic 


Pressure regulator 1 Relief Value 
Adjys^able Set to 200 psig 


Precool Line 



Gas pressurizer 
initially 2500 psig 
Decrease during usage 

Dark blue lines: liquid nitrogen precool lines 
Red lines: nitrogen gas lines 
Green: mixed venting gas lines 
Blue lines: liquid nitrogen transport lines 
Black lines: air lines 


Gas flow rate meter 


Heat exchanger 
Electric heating 


Pressure control valve 
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Test Section 

• 2 stations, 1 PT and 3 TC’s placed at each station 

• At each station a TC is placed on the top, the side, and 
the bottom of the tube 

Vacuum Chamber 



PT = pressure transducer 
TC = thermocouple 




Shaping the Future of Aerospace 


Film Boiling Correlation 


Nu fl? = h p D / k v = (7.55 x 10 4 - 7.43 x 1(T 6 z / D)Re v 094 ' (l-x e )“ 5 ' 23 Pr t 04 + 0.0568 (k, / k v )We : 0 3 Jb 

• Re^GDI/jy accounts for the single-phase vapor convection along the 
wall 


• (7.55x iO‘ 4 -7.43x I0 _6 z/D) accounts for the axial dependence of 
the convection as a cause of developing flow near the inlet 

• (1-x e ) accounts for the size of the vapor layer and thus the actual 
speed of the vapor layer. 


Pr v =(C p pi/k) v accounts for the differences between fluids. 

(k|/k v )We z 0 fb 3 , where We = G 2 zl{p,o) and 0 fb = (300 - 7J/(300 - T wet ) 
accounts for the heat flux enhancement due to liquid droplet 
impingement 
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Film Boiling Correlation 


• Model compared 
with 21,907 data 
points of vertical 
downflow data 

• MAE = 47% 
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Nucleate Boiling Correlation 


The Chen* Correlation is based on the physics of the heat transfer process and thus doesn’t involve an 
empirical correlation in terms of the Boiling number Bo = q”/(Gh fg ) 


HTC is the summation of forced liquid convection heat 
transfer and nucleate boiling heat transfer 

Forced convection with the liquid 
Heat transfer from nucleate boiling 


K= h f c +h b 


hf C ~ 0.023[/te, (l-x e )]° 8 


k =0.00122 


i 0 . 79 ^ 0.45 0.49 

K l p,l Pi 

0.5 , 0.29 „ 0.24 „ 

Pi y iv Pv 


a»[T„~T„(P)l [P,AKhP] ’S 


Accounts for the flow structure effect on convection - 
function of quality and ratio of vapor/liquid properties 


F = 


n 

U 


x 0.736 

+ 0.213 

J 


Accounts for the change in the vapor bubble 
temperature due to the thermal boundary layer 


5 = 


1 

1 + 2.53x10 ~ 6 ReP 1 F lA625 


Martinelli parameter X ff =[(l -x e )/ x e J 9 (p v / p,) 05 (^ / ju v ) QA 

14 

*Chen, J. C., “Correlation for Boiling Heat T ransferto Saturated Fluids in Flow,” Industrial and Engineering Chemistry Process Design and Development, 
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Nucleate Boiling Correlation 

• When fitting against the data it was found that the Chen 
model drastically overestimated the heat flux. 

• This had negligible affect on the results because only a 
small fraction of time was spent in the nucleate boiling 
regime 

• For larger tube thicknesses the error in chilldown time 
would be more noticeable 

• Future work will include developing a more accurate 
nucleate boiling correlation 
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Transition Boiling Correlation 

• Transition boiling = film boiling + nucleate boiling 

• The TB correlation could use either or both of these 
correlations with an additional empirical constant. 

• q” Nuc » q’Viim, so a fit was applied to the data with the 
Chen correlation HTC and the nondimensional 
temperature 6^ - (T we f— T w )l( T we f — T sa l) . 



0 . 523 er\ b 
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Rewetting Temperature Correlation 


The De Salve and Panella* correlation was used to 
calculate T wet 


T„, =T„,(P) + 0.29 d ' [T m - T m (/>)] (l + 0.2790“’ ) 


d = exp 3.06 x 10 6 / (kpC p j 


w 


erfc 


1751.5 (kpC p ) 


0.5 


w 


T m = 0.8447’ 


cr 


MAE = 5.08% 


*De Salve, M.and Panella B., “Analytical Model for Bottom Reflooding Thermal-Hydraulics in Circular Ducts and Comparison with Experimental? 
Results,” International Centre for Heat and Mass Transfer Seminar on Nuclear Reactor Safety Heat Transfer, Hemisphere Publishing Corporation, 
Managua, Nicaragua, CA edited by S. G. Bankoff and N. H. Afgan, 1982, pp. 742-762. 
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Critical Heat Flux (CHF) Correlation 


• A similar form to the Katto* correlation was used to 
calculate CHF 


q CHF = 0.0527 Gh lv We~ 02m 




per 


• MAE = 12.6% 

*Katto, Y. and Kurata, C., “Critical Heat Flux of Saturated Convective Boiling on Uniformly Heated Plates in 1 ! 
Parallel Flow International Journal of Multiphase Flow, Vol. 6, 1980, pp. 575-582 
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Temperature at Onset of Nucleate Boiling (T 0 nb) 


• There is a minimum superheat (T w - T sat ) > 0 required 
to achieve fully-developed nucleate boiling 

• If T w < T 0NB the heat flux is mostly from single-phase 
liquid convection with the w, even though T w > T sat . 

• A simple empirical fit to the data was used to correlate 
the ON B temperature: 

T 0 m = T sat (P) + 0.007 IP + 5 

• Future work will include a more robust, non- 

dimensional correlation for T 0NB 19 
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Simulation Logic 
For each solid node: 

1 . Calculate T wet for the mass flow rate and pressure. 

2. If T w > T wet then use film boiling correlation 

3. If T w < T wet . then: 

1. Calculate T onb 

2. If T w > T 0 nb. then: 

1. Calculate CHF 

2. Calculate q nb ” prediction from Chen correlation 

3 . If q nb ” > q CHF ” then use h tb (transition boiling) 

4- If q nb ” < q C HF" then use h nb 

4. If T w < T onb then use then use single-phase liquid Dittus-Boelter heat 
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1 D Simulation Description 


Only models the energy equation of the pipe material. 

d T d 2 T \ ■ / \ 

Pw C p,w = K ~rr + Iconv \ Z ) + 9 parasitic ( Z ) 


dt 


dz 


• Assume constant mass flow rate and time-independent pressure 
distribution of the fluid 

■ therefore the momentum and continuity equations do not need to 
be solved. 

• Radiation and gas conduction are included in the term q” parasit i C 

• Time integration: 1 st order fully implicit 

• Spatial discretization: 2 nd order central difference 

• L = 22.5”, N = 40, so dz=L/N=0.5625” 
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1 D Simulation Description 


At each time step the 
enthalpy of a fluid 
node, h h was 
calculated by 


^ = /v, + q 


conv A -1 
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Convergence Study 
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Simulation vs. Data 


• The simulation was run for all 54 tests 

• The usefulness of the correlations was determined by 
the accuracy of the chilldown time. 

• Two different time errors were considered 

■ the time from start to T wet 

■ the time from start to T 0NB 
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Simulation vs Data 

G = 1 26 kg/m 2 s, Re, = 1 1 ,046, P jnlet = 1 76 kPa. 
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Simulation vs Data 

G = 220 kg/m 2 s, Re, = 21 ,695, P jnlet = 252 kPa. 
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Simulation vs Data 

G = 342 kg/m 2 s, Re, = 37,810, P mlet = 420 kPa. 
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Simulation vs Data 

G = 627 kg/m 2 s, Re, = 65,575, P jnlet = 561 kPa 
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Simulation vs Data 

• G = 888 kg/m 2 s, Re, = 91 ,433, P inlet = 690 kPa. 
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Simulation vs Data 

G = 1 179 kg/m 2 s, Re, = 1 13,762, P mlet = 723 kPa 
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Conclusion 


• The ID simulation was able to predict the transient 
temperature of the tube with good agreement with the 
data over a wide range of conditions. 

• This supports the validity of: 

■ Film Boiling HTC 

■ Rewetting temperature 
- CHF 

■ the logic to choose which correlation to use 
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Future Work 

• Develop an accurate nucleate boiling correlation and 
corresponding transition boiling correlation 

• Develop a robust T 0NB correlation 

• Extend the correlations to different flow angles 
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